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Abstract 

^ ■ Recent progress in quantum field theory and quantum gravity relies on mixed boundary 

conditions involving both normal and tangential derivatives of the quantized field. In par- 
ticular, the occurrence of tangential derivatives in the boundary operator makes it possible 
to build a large number of new local invariants. The integration of linear combinations of 
such invariants of the orthogonal group yields the boundary contribution to the asymptotic 
expansion of the integrated heat-kernel. This can be used, in turn, to study the one-loop 
^ ■ semiclassical approximation. The coefficients of linear combination are now being computed 

for the first time. They are universal functions, in that are functions of position on the 
boundary not affected by conformal rescalings of the background metric, invariant in form 
and independent of the dimension of the background Riemannian manifold. In Euclidean 
quantum gravity, the problem arises of studying infinitely many universal functions. 

0^ . In classical and quantum field theory, as well as in the current attempts to develop a quan- 

tum theory of the universe and of gravitational interactions, it remains very useful to describe 
physical phenomena in terms of differential equations for the variables of the theory, supple- 
mented by boundary conditions for the solutions of such equations. For example, the problems 
of electrostatics, the analysis of waveguides, the theory of vibrating membranes, the Casimir 
effect, van der Waals forces, and the problem of how the universe could evolve from an initial 
state, all need a careful assignment of boundary conditions. In the latter case, if one follows a 
path-integral approach, one faces two formidable tasks: (i) the specification of the geometries 
occurring in the "sum over histories" and matching the assigned boundary data; (ii) the choice 
of boundary conditions on metric perturbations which may lead to the evaluation of the one-loop 
semiclassical approximation. 

Indeed, while the full path integral for quantum gravity is a fascinating idea but remains 
a formal tool, the one-loop calculation may be put on solid ground, and appears particularly 
interesting because it yields the first quantum corrections to the underlying classical theory 
(although it is well known that quantum gravity based on Einstein's theory is not perturbatively 
renormalizable) . Within this framework, it is of crucial importance to evaluate the one-loop 
divergences of the theory under consideration. Moreover, the task of the theoretical physicist 
is to understand the deeper general structure of such divergences. For this purpose, one has to 
pay attention to all geometric invariants of the problem, in a way made clear by a branch of 
mathematics known as invariance theory [1]. The key geometric elements of our problem are 
hence as follows. 

A Riemannian geometry (M, g) is given, where the manifold M is compact and has a bound- 
ary dM with induced metric 7, and the metrics g and 7 are positive-definite. A vector bundle 
over M, say V, is given, and one has also to consider a vector bundle V over dM. An operator 
of Laplace type, say P, is a second-order elliptic operator with leading symbol given by the 
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metric. Thus, one deals with a map from the space of smooth sections of V onto itself, 



P : C°°{V, M) -► C°°(V, M), (1) 
which can be expressed in the form 

P=-g ab V v a VX -E, (2) 

where is the connection on V, and E is an endomorphism of V: E G End(V). Moreover, 
the boundary operator is a map 

B:C°°(V,M)^C°°(V,dM), (3) 

and contains all the informations on the boundary conditions of the problem. Since we are inter- 
ested in a generalization of Robin boundary conditions [2-9], we consider a boundary operator 
of the form (the operation of restriction to the boundary being implicitly understood) 



B = V N + \ 



+ S. (4) 



With our notation, Vat is the normal derivative operator Vjv = N a V a (N a being the inward- 
pointing normal to dM), S is an endomorphism of the vector bundle V, T l are endomorphism- 
valued vector fields on dM, and Vj denotes tangential covariant differentiation with respect 
to the connection induced on dM. More precisely, when sections of bundles built from V are 
involved, Vj means 

V^> 1 + 1 V, 

where denotes the Levi-Civita connection of the boundary of M. Hereafter, we assume 

that 1 + T 2 > 0, to ensure strong ellipticity of the boundary- value problem [9]. 

The case of mixed boundary conditions corresponds to the possibility of splitting the bundle 
V, in a neighbourhood of dM, as the direct sum of two bundles, say V\ and Vi, for each of which a 
boundary operator of the Dirichlet or (generalized) Robin type is also given. The former involves 
a projection operator, say II, while the latter may also involve the complementary projector, 
1 — II, and the metric of V, say H [5]: 

Bi = n, (5) 
b 2 = (i - n) [hv n + x - (r l Vi + v,r) + s] . (6) 

We can now come back to our original problem, where only the boundary operator (4) occurs, 
and investigate its effect on heat-kernel asymptotics [8,9]. Indeed, given the heat equation for 
the operator P, its kernel, i.e. the heat kernel, is, by definition, a solution for t > of the 
equation 

(^ + P)u{x,x';t)=0, (7) 

jointly with the initial condition 

lim / dx J detg(x') U(x, x ; t)p(x') = p(x), (8) 
and the boundary condition 

[BU(x,x';t)] dM = 0. (9) 
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The fibre trace of the heat-kernel diagonal, i.e. TvU(x, x;t), admits an asymptotic expansion 
which describes the local asymptotics, and involves interior invariants and boundary invariants. 
Interior invariants are built universally and polynomially from the metric, the Riemann curvature 
R a bcd of M, the bundle curvature, say Q a b, the endomorphism E, and their covariant derivatives. 
By virtue of Weyl's work on the invariants of the orthogonal group, these polynomials can be 
found by using only tensor products and contraction of tensor arguments [7,10]. The asymptotic 
expansion of the integral 

/ dx^detg TrU(x, x; t) = Tr L 2 f e~ tp ) , (10) 
Jm v ' 

yields instead the global asymptotics. For our purposes, it is more convenient to weight e~ tp 
with a smooth function / £ C°°(M), and then consider the asymptotic expansion 

„ oo 

Tr L2 (fe- tp )^ dx^Mj f(x)TTU(x,x-,t)~(4Trtr m / 2 Y,t l/2 Ai /2 (.f,P,B). (11) 
V ' Jm 1=0 

Following Ref. [8], m is the dimension of M, and the coefficient Ai/ 2 (f, P, £>) consists of an inte- 
rior part, say Ci/ 2 (f, P), and a boundary part, say Bi/ 2 (f,P,B). The interior part vanishes for 
all odd values of I, whereas the boundary part only vanishes if I = 0. The interior part is obtained 
by integrating over M the linear combination of local invariants of the appropriate dimension 
mentioned above, where the coefficients of the linear combination are universal constants, inde- 
pendent of m. Moreover, the boundary part is obtained upon integration over dM of another 
linear combination of local invariants. In that case, however, the structure group is 0(m — 1) 
[10], and the coefficients of linear combination are universal functions [8], independent of m, 
unaffected by conformal rescalings of g, and invariant in form (i.e. they are functions of position 
on the boundary, whose form is independent of the boundary being curved or totally geodesic). 
It is thus clear that the general form of the A[ j 2 coefficient is a well posed problem in invariance 
theory, where one has to take all possible local invariants built from /, R a bcd ,Q a b, K{j, E, S,T l 
and their covariant derivatives (hereafter, Kij is the extrinsic-curvature tensor of the boundary), 
eventually integrating their linear combinations over M and dM. For example, in the boundary 
part Bij 2 (f, P, B), the local invariants integrated over dM are of dimension / — 1 in tensors of the 
same dimension of the second fundamental form of the boundary, for all / > 1 [1,10]. The uni- 
versal functions associated to all such invariants can be found by using the conformal- variation 
method described in Refs. [1,7,8,10], jointly with the analysis of simple examples and particular 
cases. 

In other words, recurrence relations of algebraic nature exist among all universal functions, 
and one can therefore use the solutions of simple problems to determine completely the remaining 
set of universal functions for a given value of the integer / in the asymptotic expansion (11). The 
detailed investigation of the coefficients A\,A 3 / 2 and A 2 when the boundary operator is given 
by Eq. (4) and all curvature terms are non- vanishing is performed in Ref. [8]. One then finds 
the result (which holds for all integer values of I > 2) 

A l/2 (f, P, B) = A l/2 (f, P, B) + Tr[a l/2 (f, R, n, K, E, T, S)] , (12) 

where Ai/ 2 (f,P,B) is formally analogous to the purely Robin case, but replacing the universal 
constants in the boundary terms with universal functions, whereas a// 2 is a linear combination 
of all local invariants of the given dimension which involve contractions with T\ Our task is now 
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to derive an algorithm for the general form of au 2 , since it helps a lot to have a formula that 
clarifies the general features of a scheme where the number of new invariants is rapidly growing. 
Indeed, from Ref. [8], we know that, in cti, only one new invariant occurs: JK i jT l Y^ whereas in 
a 3 /2 11 new invariants occur, obtained by contraction of P with terms like (tensor indices are 
here omitted for simplicity) 

fK 2 , fKS, fVK, fVS, fR, fSl, f. N K. 

In a 2 , the number of new invariants is 68: 57 involve contractions of with terms like 

fK 3 , fK 2 S, fKS 2 , fRK, fQK, fEK, fRS, fQS, fKVK, fSVK, 

fKVS, fSVS, fVVK, fVVS, fVR, /Vfi, fVE, 
10 local invariants involve contractions of T l with contributions like 

f, N K 2 , f, N KS, f, N VK, f. N VS, f, N R, f. N n, 

and the last invariant is f-NN^jF 1 ^ [8]. It is thus clear that the knowledge of all local invariants 
in ai/2 plays a role in the form of a^ +i y 2 , and one can write the formulae 

a 1 = / i>p)/f), (13) 
i=i 

as/2 = f±Ur' 3/2) 4 3/2) + /;ivI>P' 1) (14) 
i=i i=i 

« 2 = ft"?* 4 2) + f,N±Ur /2) 4 3/2) + f;N N ±U?» I?. (15) 

i=l i=l i=l 

With our notation, i\ = 1, %2 = 10, i% = 57, and IA[ X,V ^ are the universal functions, where i is an 
integer > 1, x is always equal to the order 1/2 of aj/2, and y is equal to the label of the invariant 

(v) 

I- , which does not contain / or derivatives of / and is of dimension 2y — 1 in K or in tensors 
of the same dimension of K. 

These remarks make it possible to write down a formula which holds for all / > 2: 

a l/2 (f,R,n,K,E,T,S) = ^f^ ^ U ^ l ' r)/2) [T 2 ]lf- r),2 [R,n, K, E, T, S], (16) 

r=0 i=l 

where f^ is the normal derivative of / of order r (with = /), and square brackets are 
used for the arguments of universal functions and local invariants, respectively. The equations 
(12) and (16) represent the desired parametrization of heat-kernel coefficients with generalized 
boundary conditions, provided that the P are covariantly constant [8]. 

One has now to evaluate the universal functions in the general formulae for A s / 2 , A 2 , A 5 / 2 
and so on. For the coefficients A^j 2 and A 2 , results of a limited nature are available in Ref. 
[8], which show that all universal functions are generated from Vl + T 2 and ^ 1 p2 ArtanhV— r 2 . 
Upon completion of this hard piece of work, one could perform the evaluation of all universal 
functions for A 5 / 2 (f, P, B) as well, possibly with the help of computers. For this purpose, one has 
to combine the conformal-variation method with the analysis of simpler cases. As shown in Refs. 



4 



[8,9], one then obtains a quicker and more elegant derivation of the coefficient A\(f, P, B). There 
are thus reasons to expect that, in the near future, all heat-kernel coefficients with generalized 
boundary conditions may be obtained via a computer algorithm in a relatively short time. This 
adds evidence in favour of the understanding of general mathematical structures being very 
helpful in providing the complete solution of difficult problems in physics and mathematics. In 
particular, from the point of view of quantum field theory in curved manifolds, this would mean 
an entirely geometric understanding of the first set of quantum corrections to the underlying 
classical theory, with the help of invariance theory [1], functorial methods [10] and computer 
programs. 

In Euclidean quantum gravity, however, if one uses the de Donder gauge-averaging func- 
tional, and if one requires invariance of the whole set of boundary conditions under infinitesimal 
diffeomorphisms on metric perturbations, one finds boundary operators of the kind (5) and (6), 
where the matrix V 2 commutes with S but not with P [8] . This implies in turn that there exist 
infinitely many different tensors of the type [11] 

a (m) (r 2 )r/3 (m) (r 2 )F'], 

which can contribute already to the integrand for A\(f, P, B), upon contraction with Kij. Thus, 
for the boundary operator given by the direct sum of Eqs. (5) and (6), even the A\ coefficient 
is unknown. 

One thus faces a highly non-trivial problem. On the one hand, analytic results exist for 
the A2 coefficient with boundary operator (5) and (6) in the particular case of a portion of flat 
Euclidean background bounded by a three-sphere [6,12]. Moreover, it has been shown in Ref. [5] 
that the boundary operator given by the direct sum of (5) and (6) leads to a symmetric operator 
on metric perturbations. However, in the non-commuting case relevant for Euclidean quantum 
gravity, even the building blocks of geometric invariants involving P are unknown. This is why 
it remains unclear how to write a general and unambiguous formula for heat-kernel coefficients. 
The solution of such a problem is of crucial importance in quantum gravity for the following 
reasons: 

(i) to improve the understanding of BRST invariant boundary conditions [13]; 

(ii) to obtain an entirely geometric description of the one-loop divergences in quantum gravity 
and quantum supergravity [6]; 

(iii) as a first step towards the quantization in arbitrary gauges on manifolds with boundary; 

(iv) to clarify the differences between Yang-Mills fields and the gravitational field; 

(v) to complete the application of the effective-action programme to perturbative quantum 
gravity. 
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